Phys. Rev. Lett., in press. 

Adsorption of Xe atoms on metal surfaces: New insights from first-principles 

calculations 
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The adsorption of rare gases on metal surfaces serve as the paradigm of weak adsorption where 
it is typically assumed that the adsorbate occupies maximally coordinated hollow sites. Density- 
functional theory calculations using the full-potential linearized augmented plane wave method for 
Xe adatoms on Mg(0001), Al(lll), Ti(0001), Cu(lll), Pd(lll), and Pt(lll), show, however, that 
Xe prefers low-coordination on-top sites in all cases. We identify the importance of polarization 
and a site-dependent Pauli repulsion in actuating the site preference and the principle nature of the 
rare-gas atom-metal surface interaction. Copyright (2002) by The American Physical Society. 
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The adsorption of rare-gas atoms on metal surfaces 
represent prototypical physical adsorption systems [1,2]. 
Noting that these adatoms have closed shells, it is typ- 
ically assumed that the adatom-substrate interaction is 
determined by an interplay of van der Waals attraction 
and Pauli repulsion, and that these adsorbates prefer 
maximally coordinated sites [1]. However, the concept of 
the van der Waals interaction (dispersion forces) is valid 
for large distances where orbitals almost do not over- 
lap, and at closer distances, near the equilibrium geom- 
etry of an adsorbate on a surface, the direct interaction 
of adsorbate and substrate orbitals can be significant. 
Thus, although the van der Waals tail is not described 
by state-of-the-art exchange-correlation functionals, the 
short-range attraction, as well as the Pauli repulsion, is 
treated properly. This appears to be valid in particular 
for the adsorption of heavier rare-gas atoms, but also for 
the scattering of He and Ne atoms at metal surfaces, the 
interatomic distances at the turning point are sufficiently 
small that density-functional theory (DFT) calculations 
provide the correct description [3]. 

With regard to the preferred adsorption site, about 
ten years ago it was argued that Xe on Pt(lll) ad- 
sorbs in low-coordination on-top sites instead of the ex- 
pected high-coordination hollow sites [4]; a quite sur- 
prising, not widely accepted, and certainly not under- 
stood suggestion. This assessment was based on He- 
atom diffraction data from the low-temperature (T < 
60 K) uniaxially compressed phase of Xe on Pt(lll) [5,6]. 
Later, a spin-polarized low-energy electron diffraction 
(SP-LEED) study of Xe-adsorption "showed" that Xe 
atoms adsorb preferentially in hollow sites on Pt(lll) [7] 
and Pd(lll) [8] in the (V3 x V3)i?30° structure (from 
now on labeled V3)- These latter findings were subse- 
quently contradicted by recent LEED studies, which re- 
ported that Xe in the V3 structure on Cu(lll), Pd(lll), 



and Pt(lll) adsorbs in the on-top sites [9]. 

Thus, despite the conceived "simplicity" of these ad- 
sorption systems, even the most basic questions have not 
been answered satisfactorily: What is the adsorption site 
and why is it what it is? and what is the nature of the 
adsorbate-surface bond? 

With respect to first-principles studies, the first good- 
quality calculation was performed by Muller [10], who 
studied Xe on Pt(lll) using the cluster approach and 
the local-density approximation (LDA) . It was concluded 
that the on-top site was energetically preferred over the 
hollow site by 30 meV. The interaction was described 
in terms of a hybridization of Xe bp orbitals with metal 
<i-states, with the on-top site affording the greater mix- 
ing and hence being the preferred adsorption site. Some 
questions about this work were raised (e.g. [9]): (£) Be- 
sides the general concern about the validity of the cluster 
approach, such calculations may only relate to low cover- 
ages, and not to the case of Xe in the y/3 adlayer, which 
is the only structure for which good experimental data 
exists, (it) It is also unclear why the top site affords the 
greater mixing. As it stands, this is more a statement 
than an explanation. (Hi) Because on-top site adsorp- 
tion was also found for Xe on Cu(lll), for which the 
(i-band lies well below the Fermi level, a crucial role of 
substrate d-states seems unlikely [9]. A recent work by 
Betancourt and Bird [11] find that Xe also prefers the 
on-top site in the V3 adlayer structure on Pt(lll). And 
for an artificial c(2 x 2)-Xe/Ag(001) structure, Clarke 
et al. [12] reported a theoretical study which also favors 
the on-top site. However, neither Betancourt and Bird, 
nor Clarke et al. presented a consistent explanation for 
the on-top site preference. Finally, a recent study of the 
Xe/Cu(lll) system using a ten- atom cluster model, has 
investigated the interface dipole and work-function, but 
the site preference was not addressed [13]. 
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1. DFT-LDA results for: (a) Relative adsorption 
energy, A Sad, of Xe on Pt(lll) (the energy of the on-top 
site is the energy zero). The inset shows a top view of the 
fcc(lll)-(l x 1) surface unit cell, indicating the six different 
sites considered for Xe-adsorption. Full circles denote the sites 
for the Xe atoms, the large open, pale, and dark gray circles 
denote the metal atoms in the first, second, and third layers, 
respectively, (b) Work-function change, A<E>. The calculated 
value for the clean surface is $ pt ( ln ) = 6.06 eV. 



To determine the adsorption site preference of Xe 
adatoms on metal surfaces, to explain the site pref- 
erence, and to provide an improved understanding of 
rare-gas adsorption on metal surfaces in general, we 
performed DFT calculations using the all-electron, full- 
potential linearized augmented plane wave (FP-LAPW) 
method [14,15] for Xe adatoms on the Mg(0001), Al(lll), 
Ti(0001), Cu(lll), Pd(lll), and Pt(lll) surfaces. We 
employed two different exchange-correlation functionals, 
namely, the LDA [16] and the generalized gradient ap- 
proximation (GGA) [17]. The surfaces are treated by 
the supercell approach using the repeated slab geome- 
tries (five- and six-layer thick slabs separated by a vac- 
uum region of 18 A). To determine the potential-energy 
surface (PES) of Xe on the metal surfaces we chose six 
sites on the fcc(lll) surface: on-top, bridge, fee, hep, 
the midpoint between on-top and fee, and the midpoint 
between on-top and hep (see inset of Fig. 1). For each 
site the height of the Xe atom and the position of the 
topmost two substrate layers are fully relaxed. 

Our calculations show: (i) The on-top site is energeti- 
cally favored for all studied systems in the -\/3 structure. 
(ii) Calculations for low Xe coverage, 0x c = 1/4 and 1/9, 
show that also in these cases Xe atoms prefer the on-top 
site (e.g., for Xe/Pd(lll) and Xe/Pt(lll) the on-top site 
is favored over the fee site by 20.59 meV and 57.87 meV, 
respectively, at 0x c = 1/9 using the LDA). Thus, lateral 
Xe-Xe interactions are not responsible for the on-top Xe 



TABLE I. Adsorption energy of Xe in the on-top site, _E a d; 
adsorption-energy difference between on-top and fee sites, E* A 
relative to the on-top site; work-function 



= El 



pon- top 



change, A"l>; vertical distance of Xe to the first metal layer, 
dxc-mctai (given as the average distance for the topmost 
layer), and substrate rumpling (a positive sign indicates an 
outward and a negative sign an inward displacement). In 
parenthesis are experimental results [1,9]. 





Sad A Sad 

(meV) (meV) 


A$ (eV) 
on-top fee 


Mg(0001) 
Al(lll) 
Ti(0001) 
Cu(lll) 
Pt(lll) 
Pd(lll) 


-130 4.34 
-176 1.12 
-196 21.82 
-277 (-190) 9.14 
-367 (-320) 49.09 
-453 (-360) 51.22 


-0.22 -0.16 
-0.38 -0.33 
-0.87 -0.70 
-0.96 -0.83 
-1.36 -0.98 
-1.44 -1.28 




dxc- metal (A) 

on-top fee 


substrate rumpling (A) 
on-top 


Mg(0001) 
Al(lll) 
Ti(0001) 
Cu(lll) 
Pt(lll) 
Pd(lll) 


3.80 3.84 
3.69 3.72 
3.56 3.66 
3.25 (3.60) 3.31 
3.07 (3.40) 3.19 
2.85 (3.07) 2.86 


0.00 
-0.01 
+0.06 

+0.03 (-0.01 + 0.02) 
+0.04 (-0.01 + 0.02) 
+0.05 (-0.02 ± 0.03) 



site preference. (Hi) Spin-orbit coupling corrections were 
included for the valence-states using the second varia- 
tional method [14]; this yields a change in the energy dif- 
ference between on-top and fee sites of < 2 meV, and the 
geometrical parameters are unchanged, (iv) The energy 
difference between the on-top and fee sites is smaller in 
the GGA than in the LDA (e.g., for Xe/Pd(lll) the on- 
top site is favorable by 13.48 meV (GGA) and 51.22 meV 
(LDA), while for Xe/Pt(lll) the on-top site is favor- 
able by 4.80 meV (GGA) and 49.09 meV (LDA)). Fur- 
thermore, the GGA predicts a weaker adsorption energy 

76 meV (GGA) and 



i t t~\ A t cXe/Pdflll) 

than LDA (e.g., E ad ' y ' = 
-453 meV (LDA), £* e / pt ( m ) = „ 8 2 meV (GGA) and 
—367 meV (LDA)). However, both functionals gives the 
same qualitative picture for the interaction between Xe 
atoms and metal surfaces. In view of these four points 
we limit our discussion in the following to the LDA re- 
sults, neglecting spin-orbit coupling corrections for the 
valence states, and discussing only the y/3 structure. All 
our results, including convergence tests, will be presented 
in detail in a forthcoming publication [18]. 

Figures la and lb show respectively the relative bind- 
ing energy (with respect to the on-top site geometry) and 
the adsorption induced work-function change for Xe on 
Pt(lll). It can be seen that the on-top site is energeti- 
cally favorable and Xe adsorption induces a decrease in 
the work- function (which is greater for the on-top site), 
in agreement with experiment. The adsorption energy, 
the relative energy difference for Xe in the on-top and 
fee sites, as well as the work-function change for all the 
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FIG. 2. Local density of states for Xe on Pt(lll) for dif- 
ferent vertical positions of Xe, Z: (a) Z — 4.53 A; (b) Z = 
3.07 A (equilibrium Xe position); (c) Z = 2.70 A. 

systems are summarized in Tab. I, where it can be seen 
that this same behavior is exhibited. For the cases where 
there are experimental values of the adsorption energy 
with which to compare (Xe/Cu, Xc/Pd, and Xe/Pt), the 
agreement between theory and experiment is quite good. 

In Tab. I we also list the important structural data. 
The DFT-LDA bond lengths are about 10% shorter than 
those obtained by LEED [9], but in view of the fact 
that the adsorption energy is weak in physisorption sys- 
tems and our theory does not include vibrational con- 
tributions, the sign and size of the deviation is as ex- 
pected. The difference in bondlcngths, however, for the 
Xe/Cu, Xe/Pd, and Xe/Pt systems, for which compari- 
son is possible, are in excellent agreement between theory 
and LEED analyses. There is one deviation between the 
calculations and LEED analyses and that is in the di- 
rection of displacement of the atom below the Xe atom. 
The experimental values lie however within the given er- 
ror margins. Particularly interesting is that the vertical 
distance between the Xe adatom and the topmost metal 
layer is smaller for the on-top site as compared to the 
fee site for all the studied cases (for Xe/Pd it is almost 
equal). This is in contrast to the calculations of Ref. [10], 
and in contrast to what might be intuitively expected. 
We will come back to this point below. 

Figures 2a, 2b, and 2c show the local density of states 
(LDOS) for Xe on Pt(lll) for three distances of Xe above 
the surface. It can be seen that there is a broadening of 
the Xe states as the Xe adatom approaches the surface, 
which leads to a partial occupation of the previously un- 
occupied Xe 6s- and 5ri-states (mainly 5(i z 2 -state), and 
a partial depopulation of the Xe 5p-states (mainly the 
Xe 5p z -state). The same behavior in the LDOS is iden- 
tified for Xe on all the studied metal surfaces. We note 
that in Ref. [19], occupation of the Xe-6s resonance was 
attributed to being the origin of the Fermi level LDOS 
which renders Xe visible on Ni(110) in scanning tunnel- 
ing microscopy experiments. This is consisent with the 
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FIG. 3. Upper panels: Difference electron density, 
n A (r) = n Xc/pt (r) - n pt (r) - n Xc (r), along the (112) plane 
for Xe on Pt(lll) in the on-top site at 4.15 A above the equi- 
librium position (left), in the on-top equilibrium site (Z=Z cq ) 
(middle), and in the fcc-hollow equilibrium site (right). Lower 
panels: Work-function change, A$, (left) and dipole moment, 
fj,, (right) versus distance of Xe from the surface. The points 
are the DFT-LDA values. 

present results. 

We now turn to consider the change in electron density 
due to adsorption, i.e., the difference between the den- 
sity of the interacting system, n Xo / pt , and the sum of the 
densities of the clean surface, n pt , and a free Xe layer, 
n Xc . We find that a polarization of the Xe atom sets in 
when it is already far from the surface, where the nega- 
tive end of the dipole points towards the surface, i.e., the 
electron density is enhanced closer to the metal surface 
(see Fig. 3, upper left). This is due to the general prop- 
erty that the potential in front of a metal is attractive for 
electrons; specifically, the exchange-correlation potential 
is more attractive on the metal side of the Xe atom than 
on the vacuum side. At this large distance where the 
Xe dipole moment starts to build up, the electrostatic 
and covalent contributions are small, if present at all. 
This behavior is consistent with the pioneering works of 
Lang [20] and Antoniewicz [21]. A corresponding po- 
larization occurs in the surface, which can be described 
as the image of the Xe dipole. The build-up of these 
dipoles happens simultaneously (and continuously with 
decreasing distance from the surface) and are related to 
the partial population and depopulation of the Xe states 
(cf. Fig. 2). In Fig. 3, the work-function change and 
dipole moment is shown as a function of distance, Z. For 
all Z, the values for the on-top site are greater and the 
variation with distance is smooth and continuous. From 
the electron density difference distributions at the equi- 
librium positions (Z=Z eq ) depletions of charge are seen at 
the Xe atom (mainly 5p z ) and at the metal atom closest 
to Xe, i.e., mainly in the Pt d z 2-states, while an increase 
occurs at the surface, between these regions of depletion. 
Similar results are found for all the other studied transi- 
tion metal surfaces. 
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FIG. 4. Left panel: Calculated binding energies of Xe on 
Pt(lll) for the on-top and fee sites as a function of dis- 
tance from the surface. The filled and open circles indicate 
DFT-LDA data. Right panel: decomposition into repulsive 
and attractive contributions (see text). Z = is the center of 
the top layer. 

The differences in the LDOS and An(r) distributions 
between on-top and hollow sites described above, show 
that the Pauli repulsion and the attractive interaction is 
not the same at each site. In order to investigate this 
aspect in more detail, we calculated the binding energy 
of Xe on the metal surfaces as a function of vertical dis- 
tance. The results for Xe/Pt are shown in Fig. 4 (left). 
At all distances above the equilibrium position, the on- 
top site has lower energy. From these results, we obtain 
Xe vibrational frequencies of E v n, = 3.79 meV (on-top) 
and 4.17 meV (fee). Reported experimental values are 
3.70 meV and 3.80 meV, [2] showing excellent agreement 
with the calculated (on-top) value. 

In Fig. 4 (right), we decomposed the binding energy 
curves into a repulsive and attractive potential using 
a standard empirical formula: E(Z) — Aexp(-BZ) — 
C(Z — Zq)~ 3 [1], where Zq is the position of the image 
plane. We note that the two right-most DFT-LDA points 
in Fig. 4 (left) were not included in the fit since DFT re- 
sults are most reliable near the minimum in the potential 
well. It can be seen quite clearly that the Pauli repul- 
sion is different for the two sites and is weakest for the 
on-top site. We also used other expressions for the repul- 
sive and attractive terms, and the conclusion was always 
the same. We furthermore decomposed the DFT-total 
energy into kinetic, Coulomb, and exchange-correlation 
energies. The results suggest that the repulsion due to 
the kinetic energy is smaller for Xe in the on-top site. 
We find the same behavior as exhibited in Fig. 4 for all 
the rare-gas atom-metal systems investigated. 

The reason why, that for Xe at the on-top site, the 
Pauli repulsion is weaker, the equilibrium distance closer 
to the surface, and the dipole moment greater, is ex- 
plained as follows: On analyzing the reactivity index of 
the clean Pt surface [22] , we find that the electronic char- 



acter at the on-top region is donor-like and at the hollow 
region it is acceptor like. Thus, for a negative charge, 
at not too close distance (or for a dipole with the nega- 
tive end pointing toward the surface), the on-top region 
can easily screen this perturbation by transferring elec- 
tron density from the on-top region (i.e., from d z 2 and 
p z orbitals) to the interstitial region. Therefore, at inter- 
mediate distances the Xe atom becomes polarized more 
strongly at the on-top geometry, because the screening 
charge at the metal can build up more efficiently. This 
is why the dipole is greater. Similarily, at close distance, 
when wavefunctions start to overlap and Pauli repulsion 
sets in, the substrate can reduce this repulsion more ef- 
ficiently at the on-top site by transferring d z i and p z 
electrons to s-states. Hence Xe can get closer to the sur- 
face. 

In summary, from DFT calculations we determined 
that Xe adatoms on the free-electron-like, early and late 
transition metals, as well as noble metals, Mg(0001), 
Al(lll), Ti(0001), Cu(lll), Pd(lll), and Pt(lll), prefer 
low coordination on-top sites. The qualitative results are 
the same for LDA and GGA functionals. We find that 
the dominant mechanisms at play in these systems are 
polarization and a site-dependent Pauli repulsion which 
is weaker for the on-top site, thus stabilizing on-top ad- 
sorption. Our recent investigations for the smaller rare- 
gas atoms, e.g., Ar and Kr on Pd(lll) exhibit the same 
behavior and basic mechanisms, and we anticipate that 
the present findings may also be generally valid for other 
rare-gas/metal systems, thus presenting a new picture 
of rare-gas adsorption, where the on-top site is the rule 
rather than the exception. 
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